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Abstract 

We consider geometric variational problems for a functional defined on a curve in three- 
dimensional space. The functional is assumed to be written in a form invariant under the group of 
Euclidean motions. We present the Euler-Lagrange equations as equilibrium equations for the in- 
ternal force and moment. Classical as well as new examples are discussed to illustrate our approach. 
This new form of the equations particularly serves to promote the study of bio- and nanofilaments. 
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INTRODUCTION 



Ever improving experimental techniques in biophysics and nanotechnology have created 
great interest in one-dimensional continuum models for such slender structures as DNA, 
p..otei,.. .anotubes a„d othe. bio- a. „anoma,„e„ts H Q Q. Q. In addit.c. such 



models continue to be used in engineerinig applications to study large statical deformations 



of one-dimensional elastic structures (e.g., cables, pipelines, textile yarns) 
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30j. Vortex 



ilaments provide another target for the application of one- dimensional continuum models 



Often these models give rise to variational problems on curves in a form that is invariant 
under Euclidean motions. The corresponding (Euler-Lagrange) equilibrium equations are 
usually derived ad hoc. To be sure, there is a general theory of Euler-Lagrange equations for 
invariant variational problems [l, 2], but it is usually expressed in abstract geometrical form 
and does not seem to be widely known in the physics and mechanics literature. Moreover, 
the equations it yields are naturally expressed as high-order ordinary differential equations 
(ODEs), which are neither necessarily convenient for further analysis or numerical solution, 
nor helpful in providing insight into the problem under consideration. 

n 

Here we show that Anderson's Euler-Lagrange equations of [1] for variational problems on 
curves can be written in the form of (first-order) balance equations for the internal forces and 
moments in the structure plus equations that can be interpreted as constitutive relations. 
We believe that this form of the equations is better suited to further analysis, in particular 
in problems of rods and filaments subjected to end loads (as, e.g., in single-molecule exper- 
iments). We demonstrate the usefulness and wide applicability of the equations by a series 
of classical as well as novel variational problems. 



EQUATIONS FOR INVARIANT VARIATIONAL PROBLEMS 

Consider the variational problem for a functional / on a smooth curve C = {r{s) G 
R^se [0,L]}: 

/ /[s,r"(s),x(s)]ds ^ extr. (1) 
Jo 

Here x{^) ^ collects possible additional functions defined on the curve. We assume 
that the scalar function / is invariant under reparametrisations of the curve C and invariant 
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under the group of Euclidean motions of M^. Then the functional can be expressed in terms 
of the Euclidean invariant properties of the curve C, i.e., its curvature and torsion. For 
such problems it is possible to write down the Euler-Lagrange equations directly in terms 
of the geometric invariants, i.e., avoiding coordinates r The following proposition (first 
briefly announced in |3|] in a slightly less general form) gives the equations in the form of 
force and moment balance equations. The result is a natural convergence of lines of work 



in mechanics, physics and mathematics that can be traced back to Sadowsky 



and Singer Capovi 
variational bicomplex 



la et al. 



16| . Langer 



and, in more abstract form, to the theory of the invariant 



Proposition. Let r{s),s G [0,L], be a sufficiently smooth regular curve in M.^ with unit 
tangent vector r'{s) = t, curvature and torsion t{s). Here and in the following the 
prime denotes differentiation with respect to arclength s. In addition, let be a smooth 
function of arclength. Then the Euler-Lagrange equations for the variational problem 

^ /(x, r, X, x', r', x', x", r", x", • • • , r^'\ X^'^) ds ^ extr (2) 



can be presented in the form of (a) balance equations for the components of the internal force 
F = {Ft, Fn, Ff,)^ and moment M = {Mt, Mn, Mb)"^ expressed in the Frenet frame {t,n,b} 
(tangent, principal normal, binormal), 

F' + X F = 0, M' + u; X M + t X F = 0, (3) 

where u = (r, 0, x)^ is the strain (Darboux) vector in the Frenet frame, (b) the 'constitutive' 
equations 

Mb = SM), Mt=£,{f) (4) 

and (c) the equations 

£^,(/) = 0, 2 = 1,2,. ..,n, (5) 
with the Euler-Lagrange operator for the variable C, defined by E^ih) = d^h — {d(^ih)' + 

Note. We adopt the notation that for any vector G M'^ the triple of components 
{yt,Vn,Vh) = {v ■ t,v ■ n,v ■ b) will be denoted by the sans-serif symbol v. Equations ([3]) 
in vectorial form read F' = 0, M' + r' x F = 0, the familiar balance equations for a one- 
dimensional elastic continuum [7|. It follows that F and M + r x F are constant vectors 
in space and that |F| and F ■ M are first integrals. 



Proof. It was proven by Anderson [])] by performing the variation of the curve that the 
Euler- Lagrange equations for x and r for the problem in Eq. ([2]) are given by 

+ - T^)8^ + 81 + 2xT8r + { ^^llz^^ + 2-8'! = 0, (6) 



2(^2 _ ^2) _ + xx"\ . ^ X' ^„ 1 



r'£:^ + 2rC - + ^ ^ 8'^ + 2—8'^ - -8'^ = 0, (7) 

\ x-^ J yc^ K 

where 7i = 7Y(/) is the Hamiltonian 

q>i>j>0 ^ ^ k=l r>i>i>0 ^ 



(8) 



q>i>j>0 k=l r>i>jl 

Equations ([5]) are nothing but the standard Euler-Lagrange equations for the functions 

Xk- We now show that Eqs. ([3]), (jlj) are simply a rearrangement of Eqs. ([6]) and (^^. Consider 
first the equation for the moment and rewrite it in component form: 

M; - xM„ = 0, (9) 

- tM, + xMt = -Ft, (10) 

Ml + rMn = Fn. (11) 

Equation ([9]) with the help of the second equation in Eq. (jlj) allows us to express the principal 
normal component as 

M„ = £;/x, (12) 
which we insert into Eqs. ( ITOj) and ( fTTl) to find 

F^ = 8'^+-8l, (13) 

X 

Fh = t8^ - x8r - (^^^ . (14) 

Next we turn to the force equation ([3]), which in component form reads 

F; - xFn = 0, (15) 

F^ - rFb + xFt = 0, (16) 

Fl + rFn = 0. (17) 
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Now, it follows directly from Eq. ([8]) that 



r'Sr 



(18) 



(here we have used that = 0, by Eq. (j^])), and if we combine Eq. (ITS]) with Eqs. (ITH) and 
( fT3]) . and integrate, we obtain 



The integration constant is fixed by the boundary conditions through the integral |F| and 
can be absorbed into the Hamiltonian Ti. This defines all the force and moment components, 
and the two equations that have not been used yet, Eqs. ( fT6l) and ( IT71) . after substitution of 
the force components from Eqs. ( fT3|) . and ( fT9|) . yield Eqs. ([6]) and (JTj). 

It is clear that the above steps can be carried out in the opposite direction, i.e., by formally 
introducing new variables Ft, Fn, Fj,, Mf, Mn, according to the above expressions one 
can write Eqs. IQ and ([7]) as a first-order system. Therefore, Eqs. ([3]) and (jl]) are equivalent 



A few remarks are in order: 

(i) Equations ([3]) and (jl]) can be thought of as arising in two steps. In the first step / is 
viewed as a function of independent variables x and r, and Eqs. (jl]) are the classical 
Euler-Lagrange equations with Mb and Mt playing the role of generalised forces. The 
order of derivatives in the operators and Sr is determined by the order of derivatives 
of X and r appearing in /. The second step then is to realise that x and r are not 
arbitrary variables, but in fact the curvature and torsion of a space curve. Eqs. (j3]), 
or equivalently Eqs. (j6]) and (jTj), are then the result of expressing the variations of x 
and T in terms of variations of the curve r. Since curvature is expressed as the second 
derivative of r and torsion as the third derivative of r, Anderson's equations involve 
derivatives up to order two in and up to order three in Sr- The balance equations 
(j3]) are a rewrite of these equations as a first-order system. The components of M 
couple the equations of step one to those of step two. 

(ii) The reason for calling Eqs. (jlj) 'constitutive' equations is that it is these equations 
that contain the physics of the problem (the balance equations (j3]) do not depend on / 
explicitly). Writing the Euler-Lagrange equations in the form of Eqs. (J3]), (jl]) and (j5]) 



Ft = H + >cE^ + tEt + const. 



(19) 



to Eqs. (jSJ) and (jTj). 



□ 
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is a way of extracting constitutive equations from the functional /. Mathematically, 
Eqs. 01]) are best viewed as equations for x and r, although they need not be resolved 
for the highest derivatives of these variables. 

(iii) Equivalents of Eqs. and ([7]) have been derived many times in the literature for 



particular applications. Examples include the isotropic Kirchhoff rod ^|, the Helfrich 



rod [8f| (corrected in [9|), piezoelectric nanobelts [10|], magnetic vortex filaments 
functionals that involve either curvature or torsion or both [sl, a functional that de- 
fends on curvature only \X2l, the Sadowsky functional for a narrow developable strip 
13], a functional that depends on x, r and their first derivatives [l^, a functional that 



involves x, r, xi aiid x'l [l^; stc. However, the explosion of terms that occurs when 
and Sr are substituted makes Eqs. (jH]) and ([7j) not particularly practical either for 
analjd^ical or numerical study (for all but the very simplest functionals /). 

(iv) It may happen that the right-hand sides in Eq. (jl]) have a simpler form in some other 
variables and accordingly we may prefer to rewrite Eqs. (jl]) (and Eqs. ([5])) in these 
new terms. Let the transformation be given by 

e = e(x, r, X, r', x', x", r", x", • • •), ^^'^ = ^ = 1, 2 

/ I I I II II II \ (i) ^ I • 1 o 

r] = ri{K,T,x,^,T ,x,>c ,r ,x ,■■■), ^ ds^' ^ = 
Then the Euler-Lagrange operators are transformed by 

^«(/) = t(-^y^ + . (20) 



, ^, . . . , 



where / is the transformed / (similar expressions hold for x as in the usual case of 
Lagrangians involving higher-order derivatives). 



EXAMPLES 



We illustrate the above theory by several examples. 



The anisotropic Kirchhoff rod 



Let the curvature xfs) and torsion t(s) define the 



centreline r(s) of the rod (up to Euclidean motions). Assuming a non-circular cross-section 
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with bending stiffnesses A and B and torsional stiffness C, we can write the elastic energy 



density as 



f{x,T,(f),(j)') = (a + 6cos20)x2 + c(r + 0')^ (22) 
where a = {A + B)/A, b = {B — A) / A, c = C /2 and is the twist angle describing the 
rotation of the local material frame with respect to the Frenet frame about the tangent 
vector t = r'. With playing the role of xi, Eqs. (jlj) and ([5]) then give, respectively, 

Mb = = 2{a + 6 cos 20) X, (23) 
Mt = drf = C{t + 0') (24) 

and 



c r' + 



+ 6x^sin20 = 0. 



(25) 



Equations ([3]), fl2^ . fl2H fl2S]l constitute a system of differential-algebraic equations (DAEs) 
that can be turned into a system of ODEs by differentiation of the algebraic equations. For 
an isotropic rod {A = B) the coefficient b vanishes and a combination of Eq. (125!) and Eq. (!24|) 
gives the first integral Mt =: c = const, which allows the system to be integrated in closed 
form. In this case the equation for the angle fully decouples from the other equations and 
the centreline of the isotropic rod can be found as a minimiser of the functional / = ax^ + cr 
with a linear torsion term [4]. On the other hand, the functional / = \A>c'^ + \Ct'^ with 
quadratic torsion was proposed to model elastic strips and polymer chains 19|, l20| . It may 
be formally obtained from Eq. fl22|) by pushing one of the bending stiffnesses, i?, to infinity 
(implying 7r/2). Rods described by this functional bend only about a single principal 
axis and therefore have their material frame locked to the Frenet frame. 

For a bundle of parallel thin rods of circular cross-section of radius R the normalised 
bending energy density may be shown to equal 21 1 

f =(}- Vl - i?2x2' 



(26) 



which provides another example of an invariant functional 34] . Ref . 221] gives extensions to 
more complicated functionals for parallel bundles to which our proposition can be applied 
to derive equilibrium equations. 

The Helfrich rod. In order to study chiral effects in polymers Helfrich proposed the 
following elastic energy density with higher-order terms included 231] 



f = f{x, 1", x') 



-X +fc3xr + — X +xrj 



(27) 



where ^2,^:3, ^22, ^4 are constant coefficients. For this functional Eqs. (jl]) become 

Mb = k2>c + 2kzXT + fc22X^ + kixr'^ - A;4x", (28) 

Mt = k^n"^ + k^x^T. (29) 

These are the nonhnear constitutive equations for the Helfrich rod (expressed in the Frenet 
frame). The second equation is algebraic and can be used to eliminate the torsion r. The 
first equation is then a differential equation for x that is to be solved in conjunction with 
the balance equations. 

The Helfrich functional has been extended to sixth order, involving the first derivative of 



24|. 



torsion and the second derivative of curvature 

A rod lying in a surface. The proposition can also be used in problems of curves with 
constraints such as the constraint for a rod to lie in a surface. If this surface constraint is 
given by the pointwise condition = g{>i,T,tp, >i',T',ilj', . . .) G M*", G M™^"*^, for certain 
m, then we consider the new functional / + A(s) • g with A(s) G a Lagrange multiplier. 

The simplest example is that of a rod in a plane. One may constrain the centreline 
r = [x, y, z) to a plane by imposing, for instance, z = 0, as in [25], but this constraint is not 
Euclidean invariant and therefore not of the type g above. A Euclidean invariant form is 
simply r = 0. We can account for this constraint by modifying the function in Eq. (l22l) and 
considering /i = / + A(s)r (hence Xi = 0; X2 = A). Equation (123|) does not change while 
Eq. fl2^ now becomes Mt = C(r + 0') + A and may be used to find the reaction A. The 
binormal force component is constant by virtue of Eq. (fT7|) . The remaining five Eqs. (ITSll . 
([16]), (El), ([To]), ([ll]) plus Eq. ([25]) with r = and x substituted from Eq. ([23]) form a system 
of six differential equations for the six variables Mt, Mn, Mt, Ft,Fn, (p. 

Note that the reaction A has the interpretation of a moment about the tangential direc- 
tion. The constraint may therefore be realised by applying a distributed twisting couple 
of the same magnitude. It may be approximated by a rod with multiple (in the limit - 
distributed) small whiskers perpendicular to the centreline (not unlike a caterpillar). If we 
imagine placing such a hairy rod between two parallel friction-free plates so that the rod 
itself would not get in touch with the plates then the normal reaction forces would give the 
required couples. 

This way of realising the constraint differs of course from the usual one corresponding to 
the z = condition, where the reactions are distributed normal forces exerted by the plates 

8 



onto the rod [25|]. Anyway, if one is only interested in the configuration of the rod then 
the reahsation of the constraint does not matter. In particular, for the isotropic rod the 
equations in both cases reduce to those of the Euler elastica, x" + = 0, corresponding 
to the functional / = x^. 

The Helfrich rod can be similarly constrained to the plane by introducing the condition 
r = 0. A more direct way to obtain the reduced functional is to delete the torsion terms 
in the right hand side of Eq. fl27|) to obtain / = ^A;2X^ + \k22x:'^ + \ki>c''^ and consider 
the problem in M?. This functional may be useful for studying polymers synthesised at 
the interface of two fiuids. It may also have application in computer vision. In this field 
the functional / = has been proposed for shape completion [26]. The Euler-Lagrange 
equation for this functional, x"' + x^x" — ^xx'^ = 0, follows directly from Anderson's Eq. ([6]) 
[l| (the equation in |26| is incorrect). Functionals that involve the torsion may be of interest 
when one deals with completion of space curves reconstructed from their planar projections. 

Rods confined to a cylinder are relevant for buckling inside tubes and for supercoiled 
filaments, and have been studiec by imposing the coordinate constraint x^ + y^ = B?, where 
R is the radius oi^ the cylinder [22] . A Euclidean invariant form of the constraint involves 



two conditions 



28|: 



g2 := >c6" — k6' + xqxcos^ ^(xq sin cos0 — r) = 0, (30) 

where >Cq^ = R and 9{s) is an unknown function that is to be found as part of the solution. 
The modified functional / + \i{s)gi + \2{s)g2 is of the required form in Eq. ([2]) (with 
Xi = Ai, X2 = A2, Xs = ^) the Euler-Lagrange equations follow from the proposition. 

Inextensible strips. An inextensible strip is a thin shell that deforms by pure bending (no 
stretching). Its surface is therefore developable and has a single non-zero principal curvature 
xi. The normalised bending energy for a rectangular strip of length L and width 2w can be 
reduced to a single integral over the strip's centreline |29| : 



I I >cl{s,t)dtds= f{x,T],T]')ds 

Jo J-w Jo 



/(X, r^, V) = (1 + V') " —, log {\^) , (31) 

where rj = r/x. In the limit w this recovers Sadowsky's functional f{x,ri) = 
2x^ (1 + r/^)^ given in {igI where Eqs. (jl]) for this case are obtained by applying the principle 



of virtual work and making use of the variation of the Frenet frame. Since the energy density 
/ depends on derivatives of the curvature only via 1], it is convenient to apply the transfor- 
mation ^ = X, r/ = r/x. Eqs. m, (EH) then yield = £ ^{f) - lS^{f) , Mt = ^S^{f). Note 
that £it{f) = dj^f and hence Mj, + riMt = d^f. These equations were first derived in 30|. 
The complexity of the centreline-reduced functional f{>i,'r],ri') makes this the first problem 
for which the invariant formulation seems to be the only way to obtain a manageable set of 



equilibrium equations. Their extension to intrinsically curved strips was considered in |31|. 

The balance equations presented here correspond to the conservation laws generated by 
the symmetry group of Euclidean motions j^. A computational procedure for deriving 
invariant Euler-Lagrange equations (analogous to Eqs. ([6]) and (JTj)) for arbitrary finite- 
dimensional transformation groups can be found in Q| . When given the balance form these 
equations may be useful for certain problems with non-Euclidean symmetry groups. An 
example is the description of world lines of relativistic particles in Minkowski space with the 
Poincare group of isometries as symmetry group 32|, l33 |. 
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